Wigner crystallization in a polarizable medium 
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We present a variational study of the 2D and 3D Wigner crystal phase of large polarons. The 
method generalizes that introduced by S. Fratini,P. Quemerais [Mod. Phys. Lett. B 12 1003 (1998)]. 
We take into account the Wigner crystal normal modes rather than a single mean frequency in the 
minimization procedure of the variational free energy. We calculate the renormalized modes of the 
crystal as well as the charge polarization correlation function and polaron radius. The solid phase 
boundaries are determined via a Lindemann criterion, suitably generalized to take into account the 
classical-to-quantum cross-over. 

In the weak electron-phonon coupling limit, the Wigner crystal parameters are renormalized by 
the electron-phonon interaction leading to a stabilization of the solid phase for low polarizability of 
the medium. Conversely, at intermediate and strong coupling, the behavior of the system depends 
strongly on the polarizability of the medium. 

For weakly polarizable media, a density crossover occurs inside the solid phase when the renor- 
malized plasma frequency approaches the phonon frequency. At low density, we have a renormalized 
polaron Wigner crystal, while at higher densities the electron-phonon interaction is weakened irre- 
spective of the bare electron-phonon coupling. 

For strongly polarizable media, the system behaves as a Lorentz lattice of dipoles. The abrupt 
softening of the internal polaronic frequency predicted by Fratini and Quemerais is observed near 
the actual melting point only at very strong coupling, leading to a possible liquid polaronic phase 
for a wider range of parameters. 

PACS numbers: 71.30.-fh Metal-insulator transitions and other electronic transitions, 71.38.-k Polarons and 
electron-phonon interactions, 71.45.-d Collective effects 



INTRODUCTION 

As it was first proposed by Wigner last century ,1] the 
long range Coulomb interaction is able to stabilize a crys- 
tal of electrons, which eventually melts upon increasing 
the density at a quantum critical point. Experiments 
done on heterostructures Q] and quantum Monte Carlo 
simulations confirm this scenario jj, |j] ■ The presence of 
impurities is known to stabilize the crystal phase in two 
dimensions |5( . Another mechanism which could help the 
stabilization of the crystal phase is the effect of a polar 
material. As a single electron moves in a polar crystal, it 
polarizes its environment creating a new quasi-particle: 
a Frohlich large polaron, with an enlarged effective mass 
[a, I3 ■ One expects then an enlargement of the Wigner 
crystal phase. 

Interesting properties of the liquid phase in polar 
doped semiconductors arise also due to the interaction 
with the polarization, as for example, the mixing be- 
tween plasmons and longitudinal optical (LO) phonons 
[g. Such a mixing can be explained by assuming a long 
range interaction between the carriers and optical lattice 
vibrations of the Frohlich type !?| . The resulting coupled 
LO-phonon-plasmon modes (CPPMs) are found in polar 
semiconductors (n-type GaP or p-type GaAs) [l^- An- 
other interesting playground for this kind of physics is 
the surface-polaron, i.e. electron close to the surface of a 
polar crystal, which have been intensively studied espe- 
cially for intermediate electron-phonon coupling a as in 



InSb where a ~ 4.5 [ill or in AgCl where a~ 3 [T^. It 
has been also observed that the gate materials (Si02,Al2 
O3 in organic thin films in field transistors are polar di- 
electrics and the interaction between the electrons and 
the surface phonons of the polar dielectric is relevant Il3l| . 
The aim of this work is to study the stabilization of 
the Wigner Crystal phase and its properties in the the 
presence of a polarizable medium. We consider a gen- 
eral model in which the key feature is the presence of 
long range interaction which arise from direct Coulomb 
interactions between electrons and from the polarizable 
medium. 

The presence of long range interactions, high polar- 
izability and low carrier density is also a common fea- 
ture of high-temperature superconductors. Of course, in 
these materials, short range interactions and lattice ef- 
fects play an important role. Nonetheless, polarons have 
been detected by optical measurements in the antifcrro- 
magnetic insulat ing phase of both superconducting and 
parent cuprates |l4lll5l ll6!|. Moreover, some evidence of 
strong electron-phonon coupling effects has been given re- 
cently in the underdoped regime |l7l |. A new interesting 
physics is introduced when studying these materials by 
the fact that the carrier concentration can be varied from 
very low to sufficiently high density. Prediction on opti- 
cal properties and more specifically the behavior of the so 
called Mid Infrared Band (MIR) by varying the doping 
has been proposed according polaronic models |18lll9ll3lj 
as well as its interpretation as charge ordering in stripes 



[20j . A similar behavior has been also found in the optical 
properties of potassium doped Barium Bismutate|21|. 

When we consider a system composed of many in- 
teracting large polarons we are faced with the problem 
of screening of both electron-electron (e-e) and electron- 
phonon (e-ph) interactions as we increase the carrier den- 
sity. A density crossover is therefore expected when the 
doping concentration is varied so that the plasma fre- 
quency approaches the optical longitudinal phonon fre- 
quency iOLO- 

At high density, phonons cannot follow the much faster 
plasma oscillations of the electron gas and therefore they 
do not contribute to the screening of the e-e interactions. 
On the other hand, the electronic density fluctuations 
screen the e-ph interaction leading to the undressing of 
the electrons from their polarization clouds. As a con- 
sequence j_the polaronic mass renormalization is hugely 
reduced |23|. In this case, the plasma frequency of the 
pure electron gas ujp = Aire^ p/m is renormalized by the 
high frequency dielectric constant £00 
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In the high density region, the self energy has been stud- 
ied by a perturbative approach for weak e-ph coupling 
{nipoi — m) in the metallic phase[23. The validity of 
this approach is ruled by the condition oiuj^Q <C ujp/eoo, 
since in this regime cop/soo is representative of electron 
density fluctuations. In this case, the electron screen- 
ing weakens the effective e-ph coupling constant and it 
is argued that the perturbative approach is suitable also 
for semiconductors which have intermediate values of the 
bare e-ph coupling in the low doping phase. The same re- 
sults have been obtained at weak and intermediate cou- 
plings by a ground state study |23|. An approach which 
is able to span the strong e-ph coupling regime has been 
presented in ref. |24| . where an untrapping transition is 
found by increasing the density via the plasmon screening 
of the e-ph interaction. There it is concluded that there 
is no polaron formation at high density, irrespectively of 
the strength of the bare e-ph coupling constant. 

At low density, the phonon energy scale (phonon- 
frequency) exceeds the electronic energy scale (plasma 
frequency). In this limit, the phonons can follow the os- 
cillations of the slower electrons and they screen the e-e 
interaction. Thus, the frequency of the electron collective 
modes is renormalized by the static dielectric constant 
Eq. Moreover, in the case of intermediate and strong 
e-ph coupling, polarons are formed 25] so that the ap- 
propriate expression for the general renormalized plasma 
frequency becomes 
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where rUpoi is the polaron mass. In the case of GaAs 
the mass renormalization due to the e-ph interaction is 



negligible, and usually eq.|(21) is used to interpret the ex- 
perimental data with tti, the band mass of the carriers, 
in place of nipoi iiQj ■ 

In ref. |26l | an approximation is developed which allows 
to study a system of many interacting large polarons in 
the intermediate/low density regime for weak and in- 
termediate coupling strengths. The phonon degrees of 
freedom are eliminated by a generalized Lee-Low-Pines 
transformation |27l | obtaining an effective pair potential 
between electrons which is non-retarded, with a short 
range attractive term and a long range Coulomb repul- 
sive term, statically screened by £o- The role of the in- 
verse polarizability parameter rj = Soo/so is evident in 
ref. |26|. In the case of r; <C 1, which is hereafter reported 
as the high polarizability regime, repulsive interaction 
and the retarded phonon-mediated attractive interaction 
are comparable leading to a softening of the energy of 
the collective modes at a finite value of the wave vector 
k, signaling a charge density wave instability. The attrac- 
tive interaction term between the electrons plays a crucial 
role also at very low density where the ground state can 
be bi-polaronic below a certain value of the polarizability 
parameter |3!tI | , or can undergo a solid/liquid phase tran- 
sition similar to the Wigner Crystallization (WC) |28l |. 
In ref. |22ll2fl|, a Large Polaron Crystal (LPC) is studied 
using a path-integral scheme. In ref 30], for ry = 1/6 
(high polarizability regime), the authors conclude that 
in the weak and intermediate e-ph coupling regimes at 
T = the LPC melts toward a polaron liquid, but in 
the strong coupling regime a phonon instability appears 
near the melting. The authors argue this behavior from 
the softening of a long wavelength collective mode due to 
the e-e dipolar-interaction. A study which shows that the 
presence of long range order is not necessary for this kind 
of scenario has been presented in ref. |^ for a simplified 
model of a classical liquid of interacting dipoles, which 
are the polarons treated a la Feynman. The dipolar 
mode (internal frequency) is renormalized by the mean 
field of the other dipoles and it is shown to soften as the 
density increases, leading to the dissociation of the dipole 
(polaron). 

The present work generalizes the approach of ref. |30j 
using a formalism which allows to span from high (77 <C 1) 
to low polarizability regime {rj ~ 1). We calculate the 
boundaries of the solid phase in three as well as in two 
dimensions. We also calculate, within the solid phase, the 
correlation function between the electron density and the 
charge polarization density. Our results confirm the rel- 
evant role of the parameter rj in the strong e-ph coupling 
regime. According to the values of this parameter two 
distinct behaviors are found: 

i) the high polarizability regime in which we found a 
scenario similar to that of ref. JSO"] i.e. the melting of 
the crystal is driven by the instability of the internal 
polaronic mode. Interestingly our more quantitative pre- 
diction push the instability-driven melting toward very 



strong couplings leaving the possibility of a liquid pola- 
ronic phase for a wider range of parameters. 

ii) the low polarizability regime, studied here even at 
strong coupling, in which we found that the undress- 
ing transition argued in the liquid phase |24| occurs also 
in the solid phase. Nonetheless e-ph interaction is able 
to stabilize the crystal against the liquid phase even for 
moderately polarizable mediums. 

This paper is organized as follows: in the first section 
we illustrate the model and the approximations used, we 
introduce the quantities of interest, and we also discuss 
the Lindemann criterion used to determine the transition 
temperature. In the second section, we present the re- 
sults in the three dimensional case. In the third section, 
the results of the two dimensional case are compared to 
the those in 3D. The conclusions are reported in the last 
section. Appendices contain technical details of the cal- 
culations. 
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I - THE MODEL AND THE METHOD 



a) The model 



The model describes a system of N interacting elec- 
trons in a _D-dimensional space, which are coupled to 
longitudinal (undispersed) optical phonons. The Haniil- 
tonian of the model is a generalization of that introduced 
by Frohlich for a single large polaron [33| to N-large po- 
larons [23| . We consider electrons as distinguishable par- 
ticles. This approximation is justified inside the solid 
phase, where the overlap between the wavefunctions of 
different localized electrons is negligible '33]. Using the 
Path Integral technique 34] phonons can be easily traced 
out taking advantage of their gaussian nature and we end 
up with the following partition function |22|: 



Z = 



jWv{n(T)\. 
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where ^ means the functional integration over all cyclic 
space-time paths of the particles fj (r) between zero and 
/3 = h/kBT. The effective electron action reads 
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Here e^ is the electron charge, m is the electron band 
mass and V is the volume, (ep,/) is the static jcllium 
charge density. The integration of phonons leads to the 
appearance of retarded e-e interaction terms -eas. (|7l8|l -. 
where the phonon propagator is 



Do(r) 



cosh{uLo[l3/2-T]) 
sinh {Pujlo/2) 



(10) 



Using polaronic units (p.u.) {hujLo for energy, I/wlo foi' 
imaginary time r and y^h/muiLo for lengths) Se-ph-e 
becomes proportional to the dimensionlcss e-ph coupling 
constant a defined as 
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while Se-e will be proportional to the e-e coupling con- 
stant 



V2e^ 
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the ratio a^/a — 2/ {\ — rj) is thus solely determined by 
r] — £oo/£o '■ when 77 ~ 1 the Coulomb repulsion over- 
whelms the attraction mediated by phonons, while they 
become comparable for 77 <C 1. Therefore, in the Frohlich 
model, the inverse polarizability parameter rules the rela- 
tive weight between the repulsive and attractive (phonon- 
mediated) interactions. This attraction can lead to a bi- 
polaronic ground state as a > adTj) |35j. Roughly speak- 
ing this condition implies strong couplings a > ac and 
high polarizability rj < rjc where ac = 9.3 and rjr = 0.131 
in 3D, ac = 4.5 and rjc — 0.158 in 2D case [35|- We 
have investigated the system for two values of rj, repre- 
sentative respectively of the high and low polarizability 
regimes, and several values of the e-ph coupling a. We 
choose respectively 77 = 1/6 as in ref.,3^], which gives 
ae/a — 2.4, and ry = 0.90519 so that the couphng ae/a 
is increased by a factor of ten [36|. For these values of r], 
no bipolaron ground state exists. 



b) The harmonic variational approximation in the 
solid phase 



We generalize the harmonic variational approach orig- 
inally introduced in ref.|29| to study the model eq.Q). 
First of all we recall here the variational theory in the 
path integral formalism. Let us consider a suitable trial 
action St which depends on some variational parame- 
ters. Substituting Seff with St in eq.Q we obtain the 
partition function Zt for the trial action and the free en- 
ergy associated to it J-'t = — fesTlnZ-r- Then the exact 
free energy can be expressed as 



where 



J^ = J^T-fci3Tln(e-*A5\ 



where AS — S^/f — St and the mean value 
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= ^fll'^mr)]{...)e--^^- (14) 



The variational free energy is obtained by a cumulant 
expansion of the logarithm appearing in eq. H13(l . At 
first order in AS it reads: 



J^v — ^T 



{AS), 



(15) 



where Tv > J' ■ To define a suitable trial action we 
proceed in two steps as in ref. [29|. First we treat the 
self interaction term S''^_ y^_^ of eq. ^ a la Feynman 



[Sllbal- Therefore we substitute Sl_f^_^ with Sp 
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(16) 
V and w are the two variational parameters. The vari- 
ational propagator Dv{t) is given by eq. H1U|) with w 
replacing lulo- We remind that Speyn ea. (|16() can be ob- 
tained by integrating out an action where each electron 
interacts elastically {Kt — rn (v'^ — w^)) with a ficti- 
tious particle of mass Mt — rn [(u^/w^) — l] . Then v 
is the internal frequency and l//i = 1/m -I- 1/Mt is the 
reduced mass of the two particle system 

As a second step, we treat the Se-e, Sj eas. (|6l9ll and 
the distinct part iS^!f^,^_J eq.® of Seff in eq.Q by 
means of a harmonic approximation. Expressing the po- 
sition of the electrons around the Wigner lattice points 
as 7\ — Ui + Xi where AT, are the vectors of the Bravais 
lattice (b.c.c. in 3D, hexagonal in 2D) and omitting the 
constant terms of the solid phase potential energy, we 
obtain the following harmonic variational action: 



St = Sk + Speyn +Sj + 6g_e + S^_^,^_^ 



(17) 
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In eg. 119(1 . the Wigner frequency is defined as usual in 
3D as w^3j5 — ujp/3 (for the 2D case see eg. ((72(1 in 
Appendix B). The force constants [XjjJ are obtained 
through a harmonic expansion for the Coulomb poten- 
tial (see appendix B). In our calculations, we neglect the 
anharmonic terms in AS of eq. I(15() , therefore we get 
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We have minimized Tv/N varying w^v sX given den- 
sity and temperature keeping a and rj fixed. Minimiza- 
tion is constrained by a convergence condition on the 
gaussian integrals appearing in Tv- The constrained 
minimization procedure is described in appendix C. 

So far, the discussed scheme appears very similar to 
the one of ref.[22|. However, we stress that the Tv, which 
we have minimized to obtain the variational parameter v 
and w , contains the hetero-interaction terms S^-e and 
which are not included in the minimization pro- 



^H,dist 
-ph — e ' 



cedure of ref. [29j. Moreover, we have also used the whole 
trial action St eg. ((171) to calculate the mean electronic 
fluctuation which we have used in the Lindemann rule, 
as explained in the following section. 



c) Lindemann rule and phase diagrams 

To determine the solid-liquid transition we use the phe- 
nomenological Lindemann criterion, suitably generalized 
to take into account the classical-to-quantum cross-over 
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in the l.h.s. of eq. 1(221) we have the Lindemann ratio 
between the mean fluctuation of the electrons around its 
equilibrium position and the nearest neighbors distance 
dji.n.- When it exceeds a critical value (r.h.s. ofeq. I(22() '). 
the solid melts. In eg. l(22() {. . .)^r, is the average taken 



over Seff eq.Q. The average is carried out at the zeroth 
order in the cumulant expansion as an average over St 
ea. ((T7|l . 

Contrary to the classical liquid-solid transition, where 
the Lindemann rule predicts the full melting line using 
a constant 7 = 7c;, in the case of a quantum crystal an 
interpolating formula for 7 is necessary to determine the 
melting line as obtained by comparing the free-energies 
of the two phases calculated using quantum simulations 
[31 . Hence the analytic expansion of the quantum correc- 
tions to the classical free energy respect to the quantum 
parameter rjq and the zero-temperature melting density 
provides the interpolating function (r.h.s. of ea l22() for 
7 iVq) |23- Vq is defined for the pure electron gas as the 
ratio between zero point and thermal activation energies 
as: 



Vq 



hcup 
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We have chosen for the function 7(77^) the form of refs. 
IS El: 



l{T,rs) = 7q 



Iq - Id 
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Formula (|24|l has a single interpolation parameter A 
which we take as A = 1.62 • 10"^ j^ 3J3 [4(i|| and A ^ 
3 • 10^2 in 2D [43. 

The chosen value of 7c; = 0.155 is such that the clas- 
sical transition lines (T = 2/rcrs a.u.) are recovered in 
both the 3D (Fc = 172 from ref. ^) and 2D {T^ = 135 
from ref. |4l|) cases. The value 7, = 0.28 is chosen to re- 
produce the zero temperature quantum transition in 3D 
{ts = 100 a.u. from ref. [40| ) and 2D (r^ = 37 a.u. from 
ref. 113). 

Roughly speaking, the transition curve is limited by 
the classical line T = (2/rc) 1/rs and the quantum melt- 
ing 1/rs = 1/r'^. The actual transition curve is a smooth 
interpolation between these two limiting behaviors. Of 
course, the precise knowledge of the interpolation for- 
mula (i.e. the knowledge of parameters appearing in it) 
is critical only for the determination of the transition line 
at high temperatures (see fig^l. 

We notice that the particular values of the parameters 
entering in eq. ^^I^ depend on the kind of statistics (bo- 
son, fermion) and on the system parameters only via the 
ratio rjq |42l | . This parameter depends on the mass of the 
particles via Wp, which measures the ze ro p oint energy 
of the oscillator which eventually melts |43|- Therefore, 
to generalize the Lindemann criterion to the interacting 
large polaron system we are left with the alternative of 
choosing between the electron and the polaron effective 
mass in eqs. (123124(1 . 

The polaron exists as a well defined quasiparticle when 
both ksT ^ hujLo |4J| and HCop ^ huj^o- The second 
condition relies on the effectiveness of the e-ph interac- 
tion, as explained in the introduction. Therefore, if both 



conditions are fulfilled, we have to replace top in eq. H23|l 
by u)p given by eq. ^. In this case, between the clas- 
sical (r/g ~ 0) and and quantum melting (r]q — > cx)), we 
have a polaronic Wigner Crystal. This is the case of the 
high polarizability (77 = 0.17). 

For low polarizability (77 = 0.9), a cross-over occurs 
inside the solid phase when hOJp ^ fiWLO f^nd the cou- 
pling is intermediate or strong, as it will be discussed 
in details later on. In this case, we still have a classi- 
cal melting of polaronic quasi-particles, but the quantum 
melting involves the undressed electrons. In the classical 
regime (low density), the transition line does not depend 
appreciably on the quantum parameter, as 7 attains its 
classical limit {i]q — > 0). In the quantum regime at high 
density and low temperatures {rjq — > 00) the function 7 
eq. (|24|) saturates to its quantum value 7gand the density 
Tc of the quantum melting does not depend of the choice 
for the quantum parameter rjq. Instead, a pronounced 
dependency on the actual value of the quantum param- 
eter is expected in the calculation of the melting line at 
high temperatures and intermediate densities. 

For 77 = 0.9 we choose the high density estimate 
lOp I ^£ca as the plasma frequency entering in eq. H23|l . 
This choice produces, in the intermediate tempera- 
ture/density region, an upward deviation (fig. ^ lower 
panel) from the classical slope. This is a drawback of our 
approximation, which is however correct at low temper- 
atures for both low and high density. 

We finally discuss to which extent we use the Linde- 
mann criterion in 2D, and more generally on the appli- 
cability of the harmonic theory in 2D. This is related to 
the well known problem of the existence of two dimen- 
sional crystalline long-range order at finite temperature 
|45j. In a pure electron gas, for T = 0, this problem 
does not arise and the properties of the system in the 
harmonic approximation have been studied extensively 
|46l l47| . The general statement for the classical impossi- 
bility of 2D crystalline long-range order was first pointed 
out by Peierls [4g. Landau [43 gave a general argu- 
ment according to which fluctuations destroy crystalline 
order possessing only a one or two dimensional periodic- 
ity. The first microscopic treatment of the problem (not 
valid in case of Coulomb interaction) is due to Mermin 
|45j : his proof is based on Bogolyubov's inequality that 
leads to the conclusions that the Fourier component of 
the mean density is zero for every vector fc in the ther- 
modynamic limit. Motivated by the interest of the 2D 
electron gas, Mermin's proof was critically re-examined 
for the long range potential [5fl, |5j| . We discuss here 
the argument of Peierls for the 2D electron crystal. The 
mean square thermal fluctuations of a generic classical 
particle diverges in two dimensions for an infinite har- 
monic crystal. At low density, we have 77^ ~ and the 
mean electronic fluctuation can be approximated by the 
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where A4_2 is the dimensionless second inverse moment 
of the density of the states (DOS) of charge fluctua- 
tion normal modes in the pure WC {p{uj)). Since long- 
wavelength acoustical vibrational modes scale as a; = c^fc, 
the DOS is given at low energies by p{uj) '^ to for 
uj ^ [53 and the integral ea. (|25|l diverges logarithmi- 
cally. However, a lower cut-off in the frequency spectrum, 
which exists for a large but finite system studied in lab- 
oratory [521 or in a computer simulation jj, |53j , removes 
the logarithmic divergence. We have chosen a cut-off fre- 
quency which corresponds to a fixed number of particles 
A^ ~ 5 • 10^. The dependence of the cut-off is discussed 
in appendix A. There and later on it is shown that our 
results are cut-off independent for low temperatures and 
density near the quantum critical point. Therefore we 
will discuss 2D case only in this region. 



d) Correlation functions and polaron radius 

We now introduce the correlation functions between 
the electron and the polarization densities for a system 
with N electrons, and a measure of the polaron radius. 

The polarization density vector of the medium is as- 
sociated to the optical phonon modes Qg- through the 
relation Q: 
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The induced charge density is defined by [g 
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Correlation between a given electron and the induced 
charge density can be defined as: 



Ci 



[r',rj 



pi{r)n^(r') 



with piif) = 5 if — fi) . In ea. H29(l we have chosen the 
appropriate normalization for the correlation function 
between one electron and the polarization. Integrating 
out the phonons we arrive at the following expression, 
in which we express all quantities in terms of averages 
weighted by the effective action eqs.Q 
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In ea. H30|l p{r' ,t) is the path density defined by: 

p{r', t)^6 (r^ - ri (r)) + ^ ,5 (P - ^(t)) (31) 

where we have explicitly separated the contribution 
pi{r',T) due to the electron 1 from the remainder. The 
first contribution in r.h.s. of eq. (|31|) give rise to a self 
term in the correlation function eq. (|30|l given by 
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Notice that in the limit of a single isolated polaron, this 
correlation function reduces to the one evaluated in ref. 
|44| . Assuming an electron at origin {f — 0), we have 
C^*^ depending only on r'. The radial induced charge 
density g{r) can be defined as 
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Using this function, we can define, as a measure of the 
polaronic radius, the square root of the second moment 
oi g{r) 



Rp — 



dr r g{r) 



1/2 



(34) 
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The actual calculation for the mean values appearing 
in eq. H32|l are carried out at the zeroth order of the vari- 
ational cumulant expansion. Explicit calculations are re- 
ported in Appendix E. 



II - RESULTS IN 3D 

Here we compare the low (77 = 0.9) and high (77 = 1/6) 
polarizability cases in 3D. For each polarizability, the 
electron-phonon coupling constant a spans from weak to 
strong coupHng regime: a = 1, 3, 5, 7, 9, 11, 13, 15. Phase 
diagrams obtained through the Lindemann criterion are 
shown in figs. ^ where the solid-liquid transition lines 
of the LPC are compared to that of the pure Wigner 
crystal. Density is expressed in term of the adimensional 
parameter rf = a^/[{A:TT/'i)p\, where ao is the Bohr ra- 
dius with (m = TOe, £00 = 1). A common feature of both 
the low and high polarizability cases is the enlargement 
in density scale of the solid phase as far as e-ph coupling 
increases. However, in both cases, the solid phase can- 
not be stabilized for any density by increasing the e-ph 
interaction, and the quantum melting point saturates at 
a maximum value when the e-ph coupling is very strong. 
To illustrate this different behavior it is worth to intro- 
duce a simplified model. 
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FIG. 1: Phase diagrams for a 3D LPC for r] = 0.17 (upper 
panel) and r\ — 0.9 (lower panel). Atomic units (a.u.) are 
used for temperature and Vs (see text). Solid phase is enclosed 
below transition lines. In both the upper and the lower panels 
continuous bold curve is the pure WC transition line and solid 
line gives the classical melting. In the upper panel dashed line 
is the renormalized classical melting. 



branch index s. The expression for the mean fluctuation 
\u^) of electrons around their equilibrium value in the 
simplified model is easily obtained by inserting w — ujlo 
in the variational expression < u^ >t (see Appendix D, 
eas ll()()ll01ll()2(l . The Sl± branches give rise to a natural 
splitting of contributions to the fluctuation 
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In the low density regime of the simplified model [30j , 
i.e. when phonons arc much faster than density fluctu- 
ations, the spectrum can be decomposed into the renor- 
malized WC frequencies J7„(cl;^ g) and the polaronic op- 
tical frequencies ri+(a; A, which can be obtained by ex- 
panding the general solutions ri|(a; ^) with respect to 
the parameter e^ j: defined as 
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(36) 



which is small for all frequencies ujk^s of WC normal 
modes at low density regime. 

The first part of the spectrum represents the low fre- 
quencies associated to the oscillation of the center of mass 
{mpoi = m + Mt) of the two-particle system, i.e. the 
electron and its relative fictitious particle (polarization), 
while the second part of the spectrum describes the dipo- 
lar modes associated to the internal motion of oscillating 
electron-fictitious particle system (fig.l of ref.|30j). Dipo- 
lar modes are weakly dispersed around the frequency ujpoi 
(eq.(25) of ref. 30J) defined as the fc = mode of the po- 
laronic branches. It represents the internal frequency of 
oscillation of the electron inside its polarization well. 



a) A simplified model 

In the simplified model, introduced in ref. [30|, the 
electrons interact with each other and with all the fic- 
titious particles ({i?i}) with mass Mt which represent 
the polarization of the medium. After integration of the 
fictitious particles, we obtain the effective electronic la- 
grangian Ceq ■ The effective harmonic lagrangian Ceq gen- 
erated by the simplified model corresponds exactly to the 
lagrangian of the action St ea. l(T7|l with the parameter 
w = lolo ■ This approximation restricts the space of vari- 
ational parameters, and therefore gives rise to a worse 
estimate for the free energy. Nonetheless, it allows to 
describe the physics of the system in a simplified fashion. 

Each WC's branch is splitted in two branches for the 
LPC and the frequencies of the system are give n by the 
two roots ri±(Wg g) (eqs. (24,25) of the work |30j), where 

uj^ j: are the WC frequencies with wave vector k and 



b) Classical and renormalized quantum melting 

Now let us consider the classical transition. This tran- 
sition is located in the low density regime of the sim- 
plified model. Using the low density expansion for the 
spectrum {Q^{uj^ ^r), il4-(a;^ ^:)), it is possible to associate 

each term, (u^), and (w^)_, to a definite degree of free- 
dom of the two-particle system, i.e. to the fluctuation of 
the center of mass 
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and to the fluctuation associated to the internal dipolar 
mode with p = u — Rt and reduced mass p 
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In this case we can easily estimate the ratio between 
electronic fluctuations in LPC and in WC by taking into 
account only the renormalized WC spectrum, i.e. the 
fluctuation associated to the center of mass eq.(j37I). Us- 
ing ea. lf^ we have < u^ >lpc / < u^ >^yc— 1/eo then 
by Lindemann criterion eq. (|22|l and by cq. H25|) at a 
given density, the critical temperature ratio also equals 
1/eo 
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(39) 



Therefore, the slope of the classical transition line is low- 
ered by the same factor, as can be seen in fig. Q] (upper 
panel), where £o is appreciably large. 

The quantum melting is ruled by the zero point fluctu- 
ations of the electronic oscillations. A zero temperature 
estimate for the pure WC gives 
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where Al_i is the dimensionless inverse moment of the 
WC DOS. If we consider only the renormalized WC spec- 
trum ea. (p?7jl . and we take into account ea. ipnji . we get 
for the LPC 
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then using Lindemann criterion we obtain at the quan- 
tum critical point (r^ = Vc) 



rc{WC) _ rupoi 



meo 



(43) 



ea. (|43|l generalizes the result of the ref. 29f| where the Lin- 
demann rule was discussed within a mean field approach. 

At high polarizability {u^) _ ea. (|37|l is the leading term 
in the mean electronic fluctuation (u^^ ea. ^Abl near 
the quantum melting for small and intermediate cou- 
plings a < 7. In this case, the quantum melting density 
scales as ea. (|^ . Notice that at weak coupling the mass 
renormalization is weak, but phonon screening through 
Eo dominates, leading to quantum melting at lower den- 
sities than in a purely electronic Wigner Crystal (upper 
panel flgs lll2|) . At low polarizability ea. (|43|l is valid up 
to a ~ 3 (fig©. 

Upon increasing the coupling, nipoi scales as ^ a"* in 
strong coupling and eq. (|15|l predicts a divergence of the 
quantum melting density. As shown in flgs. I1I2I the 
quantum melting density saturates to an a-independent 
value at strong coupling, and the prediction of eq. H43|l 
is no longer valid. We will see in the next subsection 
that deviation from the prediction of eq. H43|l arise from 
different reasons in low and high polarizability cases. 
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FIG. 2: Zero temperature phase diagram in the 2D (open 
symbols) and 3D (solid symbols) cases. In 2D q has been 
scaled according the zero density limit. Circles are the scaled 
quantum melting r^ vs e-ph coupling constant a. Dashed line 
is the renormalized quantum melting transition curve from 
ea. (l43ll . Upper panel: 77 — 0.17. Triangles locates the soften- 
ing of ujpoi- Lower panel: 77 — 0.9. The shaded area encloses 
the cross-over region inside the solid phase. 



c) High polarizability: softening of internal mode 

This is the case in which the polarization gives a large 
contribution to the total interaction energy of the system. 
The system can be thought as being composed by inter- 
acting dipoles which are made by electrons surrounded 
by their polarization. 

In the case of strong e-ph coupling we observe a satura- 
tion of the critical quantum melting density. In flg. Othe 
electronic fluctuation is reported for a — 13. Contrary 
to the small/intermediate coupling case, (u^) eQ. l|55|l 
is now the leading term near the quantum melting. The 
melting density given by ea. (|15|) is not a good estimate 
due to the contribution of polaronic optical modes which 
is now important at the quantum melting. The same 
scenario of ref.[30j is recovered: the optical polaronic fre- 
quencies drive the melting at strong coupling and high 
polarizability. Moreover we notice that (u^) ~ (l/cupoi) 
and, as density approaches a critical value, LUpoi softens 
inducing an abrupt increase of electron fluctuation which 
is dominated by the term (u^), (see figlHl- Same behav- 
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FIG. 3: The Lindemann ratio (solid line) and the function 
7 (dashed line) for rj = 0.17 and a = 13, T = 1.10"^ (a.u.). 
Contributions 5^ =< u^ >± /d^.„. of the simplified model 
ea. l35t are also shown. The inset shows the abrupt slope 
increase of the term 5i . 



ior for ujpoi is reported in ref.[30j and explained in term 
of the attractive interaction between the polarons (polar- 
ization catastrophe). We stress however that employing 
a more quantitative Lindemann criterion together with 
a self-consistent variational calculation of all Feynmans' 
parameters we get quantum melting in a region in which 
ujpoi do not actually soften. As a result the softening of 
internal polaronic frequency approaches quantum melt- 
ing only asymptotically for very large a (figEl-- 

Saturation occurs to value of r^ which seems to lie in 
the high density regime where our approach could be 
questionable. We must stress however that in the pure 
electron gas, the parameter r^ is a measure of both cou- 
pling and density. Indeed r^ can be obtained from the 
ratio of the Fermi energy to the mean Coulomb interac- 
tion, even if scaled with the band mass and the dielectric 
constant of the host medium, which are anyway fixed, i.e. 
no-density dependent. If we introduce another coupling 
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TABLE I: The critical value at the melting of the density 
(re) and coupling (r*) parameters as function of a for high 
polarizability r] = 0.16. 
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in the system, as the e-ph interaction, the two concepts 
are distinct. Global interaction is not only a function of 
the density but it is also a function of the e-ph coupling 
a. Now in the high polarizability case polarons are well 
defined as quasi-particles and we can use rfipoiipi) as ef- 
fective mass while the repulsive interactions is screened 
by £q in the low density regime. Only in this case we 
can introduce a measure of the coupling through the pa- 
rameter r* = (rnpoi/eom)rs. For the low polarizability 
case the last assumption is not valid as explained on- 
ward. The values of r* at the quantum melting (r* ) are 
reported in tab^ When the coupling a < 7 the quan- 
tum melting can be extimated thru ea. H43|) which means 
r* ~ rc{WC) = 100 that is the coupling parameter r* 
tends to the value of the Wigner crystal melting of a 3D 
electron-gas. On the contrary in the strong e-ph cou- 
pling the values of the effective coupling parameter r* 
are much bigger than those of the density Tc due to the 
huge enhancement of the polaron mass. 

Of course the exchange effects at the crystal melting 
are relevant and they can be taken into account only 
phenomenologically in our harmonic approximation (see 
footnote '42| and discussion in sec. II. C). However in the 
solid phase we must notice that these effects are ruled in 
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LPC by the parameter r* rather than r^ making them 
much more neghgible that those at the same density in 
the pure electron gas. To reahzc this fact we assume 
that the locahzed electronic wave function is a gaus- 
sian of variance a then the overlap between two of these 
functions at distance r^ is proportional to exp(— rg/4CT^). 
Now (7 in the harmonic approximation can be extimated 



as a^ = l/2mpoi(jJw where co^y = Wp^^/S is the LPC 



immediate to see that rj/ia 



1/2 a results that can 



be compared with the same for electron-gas 33] in which 
appears Vs and a different coefhcient due to a more elab- 
orate variational procedure. Taking into account data of 
table m we have that exchange effects are a fortiori neg- 
ligible in a first approximation in the case of the strong 
e-ph coupling where the quantum melting occurs at huge 
coupling parameter r* . 

In fig. 21 we show the behavior of the polaron radius as 
a function of density. While in the solid phase it remains 
almost constant, when approaching the melting density 
it suddenly increases. This behavior can be understood 
by taking into account that the polaron radius is essen- 
tially determined by the diffusion in imaginary time of 
the electron path defined in ea. H92|) of Appendix C (see 
also eas ll09ll06l in Appendix E). Its maximum value oc- 
curs at r = /3/2 which diverges at the softening of the 
polaronic frequency {copoi '^ 0). Polaronic clouds tends 
to overlap (fig0] lower panel). However, the polaronic 
nature of each particle of the LPC is preserved up to 
quantum melting. 



d) Low polarizability: cross-over in solid phase 

In this regime (ry ~ 1), the repulsive interactions 
among electrons overwhelm the attractive interactions 
due to the polarizability of the background, as can be 
seen by the relative weight of e-e and e-ph interaction 
couphng constant eas. Hllll2|) . However, self-trapping ef- 
fects are still present at least at strong coupling and at 
low density, where electrons are localized. 

Ea. H43(l quantitatively describes quantum melting in 
the low polarizability case only at weak coupling (a < 3). 
When a exceeds this value, a cross-over between a pola- 
ronic and a non polaronic phase is found inside the solid 
phase and the estimate of eq. (|43() no longer describes 
quantum melting. 

The low density regime, introduced in the previous 
subsection, is found only for the classical part of the crys- 
tal phase, where the polarization follows adiabatically the 
electron and the solid phase is a Wigner crystal made of 
polarons with an effective mass determined by the e-ph 
interaction, in the way discussed for high polarizability. 

As far as the density increases (inside the solid phase) , 
we observe that the two energy scale, lolo of phonons 
and iOp^L and eq.® of the renormalized WC frequen- 
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Wigner frequency and Wp^ is given by eq. |(5J). Then is G 
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FIG. 5: Frequencies of the system in the simplified model 
as a function of k along the direction (100), for a — 5 and 
T] = 0.9, at 1/rs = 2.10"^ and at T = 1.8 lO"'^ (a.u.). Density 
is close to the quantum melting. {Q-[ujacu],^+[u}acu]) result 
from the splitting of ujacuis,k) the acoustical WC branch. 
(f2_[LJopt], f^+ii^opt]) result from the splitting of the high fre- 
quency WC optical branch fegs. 14414511 ). For comparison are 
shown the pure WC frequencies (dotted lines). 



cies, come close and we found a cross-over region inside 
the solid phase, where the electrons and the polarization 
modes are mixed as in the liquid phase CPPM (fig.l of 
ref.|9| and fig.l of ref.[lO|). Example of the general situ- 
ation is given in fig|Sl 

To estimate the density dependence of the LPC fre- 
quencies fl±{ujs^k)j let us substitute a; r by the plasma 
frequency ujp. Results are reported in figEl which il- 
lustrates the density cross-over. In the low density limit 
(r^ -^ cx)) 17_ = wp.L eq.® while fl+ converges to ujpoi — 
V, the internal frequency for an single polaron. In this 
case, the electrons are far apart, and the " external "har- 
monic field generated by the surrounding electrons of the 
crystalline array is weak {Ke ~ e^/r^). Therefore the fre- 
quency of electron oscillation {ujp ^ Ke/m) can be lower 
than that of the phonon (lulo)j and the polarization fol- 
lows the electron oscillation. The polaron vibrates as a 
whole with a lower frequency K^/rapd ^ uip {m/mpoiSo). 
The polarization charge distribution is undisturbed as a 
first approximation, so that the value of the internal pola- 
ronic frequency {'^ v) of an electron inside its polarization 
well doesn't change . 

By increasing the density, we approach the opposite 
limit of strong external field. Now the frequency associ- 
ated to this field is too large and the polarization cannot 
follow the electron oscillation, so that each electron be- 
comes undressed from its polarization cloud. In this case 
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FIG. 6: Filled points are the typical frequencies of the sim- 
plified model obtained with lOs^h ~ ^p for a = 5 and t) — 0.9. 
Solid line is the low density renormalized plasma frequency 
eq.(|21l. Dashed line the high density renormalized plasma fre- 
quency eq.(0. Arrows mark the crossover region (see text). 



ujto is the characteristic renormalized frequency of the 
polarization. We notice that at low density r2_ gives a 
measure of the frequency of carrier density fluctuations, 
while in the opposite limit of high density, the same role 
is played by 17+. 

As we can see from figEl the cross-over amplitude is 
determined by the conditions lop,h — ^to and ujp^l — v. 
It is interesting to compare our fig|H|with the figure 1 of 
ref. Q. We notice that the asymptotic boundary given 
there by the phonon frequency ulo here it is played by 
the internal frequency. 

The cross-over of the renormalization of the plasma fre- 
quency from low to high density regime doesn't imply the 
melting of the crystal. Indeed, it is observed within the 
boundary of the solid phase estimated by Lindemann cri- 
terion. This behavior is even more clear once we consider 
the fluctuations of the position of the electrons which en- 
ter in the Lindemann criterion. 

The leading term for the Lindemann ratio at the clas- 
sical melting is S^ = (\u\ \ /d^ „ which is associated to 

the fluctuation of the center of mass ea. H37|) . Of course, 
in the classical region quantum fluctuations are ineffec- 
tive, the electrons and its polarization cloud behave as 
a single classical particle with mass rupoi. The term 

SI = {\u\')j€... associated to the internal polaronic 

frequencies ea. (|38|) is indeed negligible. 

To analyze the high density region where we meet even- 
tually the Lindemann criterion for quantum melting, we 
notice that the condition e^ j^^ 1, where e^ ^ is defined in 



eq. (|3(j|) . can be fulfilled by the majority of normal modes 
at high density. Of course, long wavelength acoustical 



and even " optical "modes in 2D have vanishing energies, 
but their spectral weight is low enough to be neglected 
in the following considerations. Expanding ^±[us^k) in 
1/e r we get 






UJLO 



(44) 
(45) 



In figOthe general solutions Q±{uJs^k) are shown for all 
the branches of the simplified model near the quantum 
melting. 

The branches {fl-[u!opt],^+['-^opt]) which results as 
splitting of optical model of the Wigner crystal ujopt{s, k) 
are well described by approximation of eas. l44l45() . 

The frequency dispersions (ri_[ci;Qcu], i^-i-[wacM]) of the 
modes which originate from the splitting of acoustical 
branches of the Wigner crystal is also reported. While 
at short wavelength, the dispersion approaches the esti- 
mates given in eas. (l44l45|l the long wavelength part of 
the spectrum is conversely described by the low density 
expansion Q,±. 

Thus we have that at the quantum melting the low 
energy part of the spectrum still behaves as in the low 
density regime. The modes depicted in the lower part of 
fig m belongs to this part of the spectrum. 

A measure of the wave- vector below which we have this 
behavior can be obtained by the condition e^ j: = 1 . The 
associated energy scale is given by uj^ — muo^Q/ {e^nipoi)- 
Contrary to low density regime eas. l|37l38|) . it is not 
possible to associate to each term of the fiuctuation 
eas. Hl()lll02|l a definite degree of freedom. However, ex- 
panding the electron fluctuation with respect to the pa- 
rameter e^ J for the frequencies uia,k < ^c and with re- 
spect to the parameter 1/e^ j: for the frequencies ujs,k > 
LOc and using eq. (|35|l the electron position fluctuations 
can be approximated by: 
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Notice that the interpretation of the fluctuations asso- 
ciated to electronic motion in this case is different from 
that valid at low density. In particular the high energy 
contribution (the second term of eq. H47|l represents a 
Wigner crystal-like fluctuation with a low energy cut- 
off. This is the largest contribution to the fluctuation at 
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quantum melting and does not depend on e-ph interac- 
tion. 

Indeed the leading term of fluctuations at quantum 
melting is (m^) , ■ This is due to the vanishing of the 
spectral weight associated to the low frequencies uj < tOc 
at high density (eq.l^Hl))- 

The saturation of the quantum melting point can be 
seen in the phase diagram of fig. Q] (lower panel). Two 
comments are needed. First, in the case of very low e- 
ph coupling, the density crossover does not occur inside 
the solid phase. Therefore, these arguments do not apply. 
The quantum melting point depends on the e-ph coupling 
as we have discussed in the previous section. However a 
saturation of the quantum melting density is observed 
clearly in fig^for intermediate and strong coupling. As 
a second point we have to emphasize that the quantum 
melting density is not that of a purely electronic Wigner 
crystal. 

This fact can be explained by writing the total electron 
fluctuation as the sum of the two terms < u^ >Hi and 

< v? >Low where < m^ >Hi is the contribution to fluctu- 
ations of modes having energies higher (lower) than tOc- 
We notice from eq (I47II that in both LPC and WC case 

< u^ >Hi are the same. But while in the WC case the 
two terms are of the same order < v? >low—< u^ >Hi, 
in the LPC case < u^ >Low^< u^ >Hi as far as the den- 
sity increases. This is due to the renormalization of the 
low energy frequencies. Therefore, the electronic fluctua- 
tion in the LPC increases more slowly with density than 
those of the WC. At given density, the electronic fluctua- 
tion of the WC is greater than those of the LPC and this 
fact explains the shifting of the quantum melting toward 
higher densities. 

The cross-over is also evident in the polaron radius. 
In the upper panel of fig. |3 we plot the polaron radius 
as defined by the eg . H34|l . We see that for any value of 
the e-ph coupling, the polaron radius tends to decrease 
as far as the density is increased. We recall that as far 
as the renormalized plasma frequency eq.|(21l exceeds the 
phonon frequency, we enter in a region in which the polar- 
ization is adiabatically slow compared to the electronic 
motion. Therefore, the electronic charge appears as a 
static distribution whose radius decreases upon increas- 
ing the density and the polaron radius follows this trend. 
The crossover is evident by scaling the polaron radius 
with Tj, as reported in the lower panel of fig. Qat inter- 
mediate and strong a. Notice that as in the high polar- 
izability case at the transition Rp/rg ~ 0.475. 

It is possible to estimate the high density limit of the 
radial distribution of the induced charge. Using the con- 
dition LOo ^ w/^Eoo we have the following expression 
valid at low temperature {ksT <C hujo) (for details see 
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FIG. 7: Polaron radius in polaronic units (upper panel) and 
polaron radius scaled with Vs (lower panel) vs (l/r^) (a.u.) for 
different a in the case rj = 0.9 at low temperature {T — 5 10~^ 
p.u.). Filled points refers to the solid phase. 
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The first term of eq. H48|l takes into account quantum 
charge fluctuations which are relevant at small distances, 
while the remaining term is a classical contribution com- 
ing from the static charge distribution. Notice that only 
the first term depends on the c-ph interaction. There- 
fore, the polaron radius tends to the same high density 
asymptotic value for different values of the e-ph coupling 
a (see upper panel of figCjI. 

As a last point we notice that the cross-over condi- 
tion, roughly estimated as ujp ^ ujlo ^ (1 — 'qf jo? ^ 
O.Ol/a^, shifts toward higher densities as the e-ph cou- 
pling constant a is reduced. In the weak coupling regime 
it lies in th e liq uid phase where RPA approaches in both 
3D refs.[2^|5J| and in 2D case [55j can be applied. It is 
also worth to remark that for high polarizability and for 
all coupling the polaronic crossover is located in the liq- 
uid phase according to the highest value of lolo ^ ^.l jo? 
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III - 2D CASE 

The results obtained in the 2D case are qualitatively 
similar to the 3D case. Both the cross-over phenomenon 
in the low polarizability case and the softening of the 
polaronic ft'equcncy in the high polarizability case are 
observed. Results are reported in the zero temperature 
phase diagram of fig. |5] In this figure, we compare the 
phase diagrams in 2D and 3D by scaling appropriately 
the 2D e-ph coupling constant following the single po- 
laron results of ref.j^: a^o = (37r/4)a2D. 2D and 3D 
melting curves scale well according to the zero density 
scaling for all studied cases. A discrepancy is found in 
the the high polarizability strong e-ph coupling softening 
of ujpoi ■ Let us first discuss the scaling at finite density. 

In our variational scheme, the DOS of the WC is the 
peculiar difference between the 2D and 3D cases. To see 
this explicitly let us compare the e-ph interaction terms 
S^'Lih-e- Assuming polaronic units we get: 



~>self 
e—ph—e 



T,3D 



3N 

oself 
e—ph—e 



Do{i 



_(a)^ /V ;^°''' - (49) 



T.2D 



2N 



37r \ V2 



Do{r) 



— a I / dr — : 

4 / 6 io V¥M^ 



(50) 



where the imaginary-time diffusion dir) ea. (j92|l is itself 
a functional of the DOS. We notice from eas. H49l50(l that 
the free energy functional scales explicitly as in the sin- 
gle polaron case [SQ by scaling the coupling constant 
a. Related to the different 2D and 3D DOS we remark 
the different behavior of the frequencies of the normal 
modes. Noticeably, the "optical" branches go to zero as 
~ Vfc at long wavelengths [^. As in the 3D case, the fre- 
quencies of the LPC are splitted in four branches (figEl) 
f^±[i^ocu(s, /c)] and n±[ujopt{s,k)] according to the same 
equation of 3D (see figlHIl, where the 2D value for lu\y is 
given in appendix B eq. H72I) . 

Let us discuss the deviation from the the scaling at 
strong coupling, which we see from fig|21in the density of 
the softening of the polaronic frequency Upoi- Actually 
we observe that a steep fall of the variational parameter 
v{rs) occurs as density increases determining the soften- 
ing of ujpoi- Peculiar features of the DOS enters in the 
variational determination of v{rs) as can achieved by the 
following argument. First of all we assume that w is 
very close to the value lulo at strong coupling. Then we 
notice that as in 3D high polarizability case the renor- 
malized plasma frequency is much less than the phonon 
frequency and the discussion which follows eq. 136|l holds 
for all densities lower than critical density of the soften- 
ing. In this case the spectrum is composed by the low 
energy branches (renormalized WC) and the by polaronic 
branches weakly dispersed around v (see also figlHJ. Us- 



ing this results at low temperatures (/3 
dition for the extrema of J^v reads: 



oo), the con- 



1 



„3/2 



— ^TWs + gia,T],r,,v) = (51) 



we, 



where the first and second terms are the derivative of 
Tt ea. H89() and g is the derivative of eas. H90l91l93() from 
Appendix C. As ti(rs) ^ for rs ~ Vc the second term ac- 
quires importance and DOS enters in the second moment 
A^2. However there are other terms which are divergent 
as w — > coming from the explicit form of the function 
9- 
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FIG. 8: 2D case. The eigenfrequencies of system along the 
direction (10) for r; = 0.9 and a = 2.12. l/r^ = 8 10"^ 
(a.u.), T = 2.5 10~^ (a.u.). Density is close to the clas- 
sical liquid-solid transition. Upper panel: The frequencies 
of polaronic branch weakly dispersed around the polaronic 
frequency ujpoi — v (dashed line). Lower panel: The renor- 
malized Wigner crystal branches (points) and the pure WC 
branches (dashed lines). 



CONCLUSIONS 

We have studied the behavior of a low density elec- 
tron gas in the presence of a polarizable medium, where 
polaronic effects may play a relevant role. To determine 
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the transition line, we have used a generahzed Linde- 
mann criterion which reproduces correctly the pure elec- 
tron gas quantum melting. Because the amplitude of 
quantum fluctuations depends on the e-ph renormalized 
plasma frequency, the Lindemann rule has been critically 
re-examined and adapted to the polaron crystal. This 
procedure allows to determine quantitatively the phase 
diagram of the model and to extend the study of the 
model to the low polarizability case, which was not stud- 
ied before. We have also studied the 2D case, showing 
that the dimensional dependence is not crucial to deter- 
mine the nature of the quantum melting within in our 
harmonic variational scheme. The scaling predicted for 
the e-ph coupling constant at zero density do apply as 
well at non zero density up to quantum melting. A no- 
ticeable difference instead is in the position of the density 
of the softening of the polaronic frequency which in 3D is 
much closer to the melting than in 2D case. This suggests 
that the hetero-interactions are less effective to destabi- 
lize the dipolar crystal in 2D. However other possible 
mechanism, (lattice-effects, structural disorder or impuri- 
ties) can cooperate to the localization together with the 
interactions between the electrons and lead to the forma- 
tion of a pinned Wigner Crystal. In this case the melting 
can not be predicted by the Lindemann rule but a sim- 
ilar dipolar instability due to the long-range interaction 
between the electrons can still drives the melting. 

While the weak e-ph coupling regime is similar for both 
low and high polarization case, the strong coupling sce- 
nario is qualitatively different. 

In the high polarizability regime, we have recovered the 
incipient instability which was found in previous studies 
near the solid phase J23, l2!j| and also in the liquid phase 
|3l| . In comparison with previous work, we have found 
that this regime is restricted to very large values of the 
coupling a > 10 leaving an interesting intermediate re- 
gion of coupling in which polarons may exist in the liquid 
phase. This region can in principle be explored with non 
perturbative numerical techniques as e.g. Path Integral 
Monte Carlo. Work along this direction is currently in 
progress. 

In the low polarizability regime, a crossover occurs in- 
side the solid phase when the renormalized plasma fre- 
quency approaches the phonon frequency. At low den- 
sity, we still have a LPC, while at higher densities the 
electron-phonon interaction is weakened irrespective of 
the bare electron-phonon coupling. In this case, polaron 
clouds overlap as and the polaron feature of the crystal 
is lost. The crossover from polaronic {cup < ujlo) to 
non polaronic {uip > lolo) character has been observed 
in weakly coupled systems such as GaAs in the liquid 
phase and analyzed in term RPA [22, |55| . In this sys- 
tem it occurs around Vs ^ 0.6 — 0.7 while for ZnO a is 
larger shifting the crossover to r^ ^ 7. Finding a sys- 
tem with low polarizability and larger e-ph coupling is 
difficult since it implies very low lu^o'- a ~ (1 — vi)/ujlo 



from ea. Hll(l . However in surfaces of InSb an a = 4.5 
has been predicted together with t] = 0.88 |11| leading to 
the possible observation of the crossover inside the solid 
phase. 

We notice also that our low polarizability scenario of 
density crossover inside the solid phase bear some resem- 
blance to that found for ripplonic polaron systems |58l |. 
Though the electron-ripplon interaction in these systems 
is different from the Frohlich type, resonances in the ab- 
sorption spectrum observed by Grimes and Adams [60j, 
their explanation at low density |6ll | relays on the same 
qualitative arguments developed in the present work. Re- 
cent works on high density ripplonic polaron systems re- 
alized on a helium bubbles predicts also in this case a 
mixing between plasmon and polaron modes |62l | . 

Finally we remark that we have obtained an apprecia- 
ble stabilization of the crystal phase even for intermediate 
regime a ~ 3 — 5 in low polarizability cases. We con- 
clude that the general result that e-ph interaction effects 
can stabilize the Wigner crystal phase could motivate ex- 
perimental studies on two dimensional electronic devices 
involving polarizable media. To this aim a layered con- 
figuration is advised even with some warnings |63| . In 2D 
heterostructure the use of a perpendicular electric field 
|63l | could not only increase the polaron effect but also 
tune it as it was shown in the case of charged helium 
surfaces l64l. 
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APPENDIX A: THE LOW ENERGY CUT OFF IN 

2D 

The 2d DOS function can be defined as 






d^k 5{u} ■ 

Vb 



UJs.k) 



(52) 



where Vb is the volume of the first Brillouin's zone 
[IBZ). Let us consider a small fraction e of the plasma 
frequency ujp. At long wavelength (fc = 0), we have the 
2D dispersion laws for the acoustical mode is wi(fc) ~ cifc 
while the "optical" uj2{k) — C2-\/fc |4g. As a consequence 
the behavior of the DOS for w ~ is 



15 



/ de I dkk 6{lo - cik) = —uj (53) 

Jo Jo Ci 

/ dO dkk 6{uj - c2\/fc) - —uj^ (54) 

Jo Jo ^2 

Introducing the scaled frequency x defined as lo = ujpx 
and the quantum parameter riq ea. (|23|l . the thermal elec- 
tronic fluctuation is expressed as an average on the DOS 
as 



/ 2\ ^ /coth(77qx) 



(55) 



DOS 



= — /dx p{ujpx) + % ( ao + a2{T]qx) + . 

(56) 
Since p{ujpx) ~ x for x ^ 0, the average in eqs. H55|) 
converges for any of x^" with n > in the the expansion 
ea. (|56|l . In the n — ~1 term we consider infrared cut-off 
Xc giving 



DOS 



dx 



{2Tilcl)ujp 



2tiujp I e 



dx 
'dx 



p{topx) 
x^ 

x^ 



(57) 



This term diverges logarithmically as Xc — > 0. How- 
ever 77g — > CXI as we approach the quantum region. The 
electronic fluctuation turns out to be cut-off independent 
if 



(58) 



^ ) < ry (ao + a2(?7gx) + 



We have chosen for the cut-off frequency Xc = 
'^min/^p — 5 • lO^'"^ so that the condition ea. H58() is 
fulfilled around riq{T,rs) > 10 which corresponds to a 
large region inside to the solid phase. By the rela- 
tion for acoustical long wave excitation ujmin = cikmin 
and kjnin — 27r/ ( r<,\/iV] , the number of electrons is 
N = 5.24 • 10®. Our inverse second moment of DOS 
is M-2 = 12.5 (cfr ref. ^] M-2 = 8.16 for N = 1024). 



APPENDIX B: THE HARMONIC VARIATIONAL 
APPROXIMATION 

We expand in the harmonic approximation the terms 
5e-e,5e-j,5*'*;^_, (eqsEEEl). Let n = R^ + u„ where 
Ri is the lattice point of the crystal and Ui is the 



electronic displacement from R^, and set Amj.,(t, cr) = 
Uj{a) — Ui{t) and _Rj.i = Rj — R^. The static terms give 



(59) 
the e-ph interaction docs not change the equilibrium 
positions of the pure electronic crystal (WC) which cor- 
responds to the minimum of the potential energy. 
The sum of the dynamical parts in the harmonic approx- 
imation gives 



S. 



H I cH I qH.dist 
e-e I '-'e-J + '-'e-ph-e 



dTy^[Vj{u,{T))+V,_,{u,{T))] 

Jo 



(60) 



where 



VAu^ir)) = -'-^ Id^r 



£0 

.2 /./3 



' -i|(61) 



Ve-eiuAr)) = ^ /daF(r-a)Vc/,,,(T,a) (62) 

Uj^T^a) = -Au,j(t, cr) •X(^j,,)Au,j(t, fj) (63) 
F{r-a) = ^-^^^^-"^D^ir-a) (64) 
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R 



■j,« 
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R 



■j,t 



(65) 



^From now we drop on the double a, r indexes in Auj j. 
To evaluate the integral in ea. H61|l and the sums on index 
J in ea. l|^ . we consider a sphere Sr of radius Rg (a disk 
in 2D) centered on site i. We first sum on index j and 
then we perform the limit Rs —> 00. Finally we sum on 
index z in cq. I|60|l . 



3D case 

By Gauss's law (with the condition Vj(0) — 0), we 
have 



^w I 



K7(^?,(r))-im^K(r)r 
z £0 



where the Wigner frequency is uj 



,2 _ 



(66) 

— - . -„ -^.^ ^ -_ - ^^ -p/S- Because 

of Vj is independent of the size of Sp, ea. l(HB|) does not 
change in the limit Rg -^ 00. 

To evaluate the sum in eg 1621 with the condition Rj < Rs 
we remind that we have two self terms {{i,i) and {j,j)) 
and two distinct terms {{i,j) and {j,i)) in eq. (|63|l . The 
two self terms give the same contribution ,as can be easily 
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check if we firstly we carry on the Hmit Rg ^ oo and then 
the sum on index j and i. They vanishes because of cubic 
symmetry of the lattice. When the two distinct terms 
{{i,j) and {j,i)) of eq. (|^ are inserted in ea. l|^ and 
the limit i?s — > cxd is taken, we obtain the term T4-e(wj) 
of eqinOI 



(67) 
Summing on index i and integrating on variable r 
the eas. (|66l67ll . we obtain the terms S^_j,S^_^,S^Zph-e 
eas. (flMn|) . 



2D case 

In 2D the interaction potential V/^(m) of a uniform 
positive charged disk of radius Rs eq. (|t)l|) is 



2 /-Sir 

VP{u) = -^^ d9F{9) 
eo Jo 



(68) 



where 



F{e) = ^Rl+u^-2RsUCoa{e) -u-R^ 

+^ cos(^) In „(l-co.(e)) 

In the limit R^ —f oo 



lim V^f (u,) = lim -pj—u^ = (69) 

since the total electric field of an infinite charged disk is 

perpendicular to the disk. 

Then we have to evaluate the sums ea. (|63f) . The two 

distinct terms ((i,j) and (j, i)) gives the identical result 

eqElof the 3D case while the self term {i,i) is written 

as 



2 



( \ 



Ui = u^'Dui 



(70) 



the matrix V in 2D is defined as sum of the matrices 
I{Rj) ea. H65l) on hexagonal lattice points Rj. Contrary 
to the 3D case, the matrix T) is not zero in 2D case. By 
the lattice symmetry, the off-diagonal elements are zero 
while the diagonal terms are equal to the local potential 
which acts on each electrons 



1 w' 



-/ daF{T - a)V^^\u,iT)\' = -m^^\u,{T)\' (71) 

2 eo 



where we use as definition 2D Wigner frequency 

^2 



Ww = — lim 
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^w 



Rj ,<R 



(2D) (72) 



For an hexagonal lattice of nearest neighbor distance 
dn.n., we have E,^»(l/2J? y = 5 51709/4.„.. 
Summing on index i eas. (|67l71|) we obtain the terms 



Normal modes 



The WC normal modes are defined as 
1 



u. 



k,s 



JkR, 



(73) 



where the vectors k belongs to the IBZ of the reciprocal 
lattice, it are eigenvector with eigenvalue uj'i of the 

^ k,s o b ^^g 

dynamical matrix A4 which is defined as 



A^Q/3 — SapUJw ^ Z_^ 1af3{Ri)e- 



ik-R, 



(74) 



-R.#0 



Inserting the WC normal modes ea. H73(l in 
eas. (|16ll8ll9l2U|l . we express the harmonic variational 
action St as 



s.k 



(75) 



where the Lagrangian is 
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mwiv —w j 
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8g 
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(76) 



APPENDIX C: THE VARIATIONAL FREE 
ENERGY Tv 

The first term of the variational free energy Tv eq. 
(|21|l is the free energy Tt associated to the partition 
function of the trial action Zt- This is calculated as the 
functional integral cq.lO where Seff eq. Q is replaced 
by St eq. (|17|) . The second term of Ty is the mean value 
eq. (|14|l of the difference between 5g^p^_g eq. (Q and 

Speyn eq. ^^■ 

We start by changing the dynamical variables of integra- 
tion from {'u.j(t)} to {g^ j:(t)}. By reality condition we 
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have a r. = ^7? and e ,- = —it i we must sum only 
k vectors in the upper half space (fcz > 0) of IBZ 
1 



N 



/ J k,s 



1k,s^ 



ikR^ * fi — ik Ri 



(77) 



s,fc,fc^>0 

Therefore the real and imaginary part of q^ j: for all k 
with (fcz > 0) of the IBZ are the actual independent 
variables and the Jacobian of canonical transformation 
is J = 2^^ 

s,k.k^>0 

(78) 

Using the periodicity condition (g^ ^(0) — q^^{(3)), 

we have the following Fourier expansion (a;„ = (27r//3)n) 

9,£„ = -p [dr q^j:{r)e-'--^ 



Is.ki^) = 


= 9s,fcx + -^^s.fe W 
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= ^i ^-M(-) 
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^'z.el-) = 




(81) 



where we have separated the mean value of path on the 
imaginary time ea. (|80|l (centroid) from the fluctuation 
around it ea. l|HTl) . The action Sxiiq^ fc(''')}) ^^ quadratic 
in {q g } therefore we can separate ea. (l78ll in two gaus- 
sian integrals 
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(r>f-^io)(»f-u;j,) 



^1 = '7o2 o2''uo2' nlV (cyclic perm. 7 =1,2, 3) 
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(87) 



the frequencies J7^ (7 = 1,2,3) are the opposite of the 
roots of cubic 

Vz(^z) = z^ + a2z'^ + aiz + ao 
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02 = w + (jJ^o H 1 : 
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The gaussian integrals eq. H84|l are convergent if A^ ^ ^ 

are positive numbers V(s, fc, n). This condition is fulfilled 
if fl^ are all positive. The numerical minimization of 
the variational free energy has been made enforcing this 
constraint. Performing the infinite product in ea. H84l) we 
have 



-■TJq 



smh{hu}Lo/'2kBT)\ f smh{hw/2kBT)^ 



n 

s,/i:,7 



hwLol'^kBT 

f^^,,srk/^kBT 
smh{hn^^^^^/2kBT) 



hw/2kBT 



-.Re rln^"^^ 



n 



TT ^k,s,c ^k,s,c i^BTZr^~77^ 

s,k,kz>0 

kBT 

r^^s,k/V^ 



(83) 



hence after we omit the classic term Zt^c ea. (|83|l . 
Zr.Sq = 



dqS" dqi"' 
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S,fc,fe2>0 



and finally we substitute the sum on (fcj, s) with the in- 
tegral on the WC DOS p{uj) in the free energy !Ft 
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To calculate the mean value of 5^1^^^ eq. iQ in 31? 
we use the following identity {57j | 
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while in 2D {q^ = q\+ql) 
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where doir) is the imaginary time diffusion in the LPC 
defined as (3D or 2D) 



(92) 
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The mean value of Speyn eq- H16() is 
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To obtain eqs. 



31 wc have used 



gig-[«,(T)-«,(CT)]\ _ g-^(iD(T-CT)g^ 



We will demonstrate ea. (|94|l in the next subsection. 

Calculation of (exp(ig • [uz{t) — Ut{a)]))j, 
^From eas. (|73l81f) we have 

iq- [U^{T) - Wj(cr)] = Y^ [<lk,s,n'^lkA'^ - CT, <?) + C 

s,fcz>0 



JlkA'' - ^,^0 = -^q- £k,s (e 



'i^ji'T ^'i^n(7'\ ^ifc-Rj 



(95) 



(96) 



then we have 

(exp(i(f • [u,(r) - Ui{a)]))j. = 

1 r dqS-" dqi"^ n. h,s.^ 
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APPENDIX D: MEAN ELECTRONIC 
FLUCTUATION 

The relation between the mean electronic fluctuation 
and the imaginary time diffusion du (r) eq. (|92|l is 

do (r) = § [< l"(0)l' > - < u(t) • u{0) >] (99) 

comparing ea. (|99|l and ea. H98|l for doir) and inserting 
the DOS function, we have 



< u > 



CTy = 



D 



= dujp{uj)Y^ 



hA^{uj) ( l3Vt^{uj) 



^^ 2m^(tj) 



coth 



(100) 
If we fix w = cjo, we have ils = ujo for one solution of 
the cubic polynomial eq. I|88(l and by eq. H87|l we have also 
A3 = 0. The other two terms give 

(u^)^^ldujp{u})-^ — 7^ TTzr^r coth ( ;^^ ) (101) 



(m^)_ ^ dojp {uj) 



nl-nl 2mf7i \2kBTj 



^'^ coth ' 



^%-^^, 2ma 



2kBT 



(102) 



Notice that if we take a single Wigner frequency being 
representative of the electronic spectrum {p{uj) = 5{lo — 
low)) we recover the results of ref. [22|. 

APPENDIX E: POLARON RADIUS 

We now calculate the density-density correlation func- 
tion of the ea. (|32|l for the variational harmonic action 
St- We assume that the equilibrium position of the ref- 
erence electron « = 1 is the origin. With the same method 
to obtain eas. (|^ . we performed the following Gaussian 
integrals for the density distribution pi (r) 



ZT,5qJ -'■-'■ nkBT/muj^ 

s,k,k^>0 



[2tt} 



[27 



.2 1^/2 



and 






k„T 



A. r Jr 



?ii:.,.|9-ej.,p.,.Jr~.)g^ 



iqui iq'{ui{T)~ui)\ _ 



^T „2 dJT) 



e 2 'J e 



Vi'-W+'Tl 



(103) 



(104) 



iTrDT,s.kdu,,jT-a)q 



^-^dD{T-a)q'^ 



(gy'jinserting ea. p04|l in ea. (|5^ . we have the density-density 

correlation function in the imaginary-time for the ^ = 1 
where the component of frequency ujs,k of the imaginary pipptron 



time diffusion duir) is [A^ = ^-^(ws^fc), 51^ — Q.^{ijJs,k)) 



s.k 



pi{r)pi{r\T))^ = (p,(f))^f__!l^:__ (105) 



{2ttP{t)) 



D/2 



1 Y^ ^ Y^ 

'nd ^m^ni 

s,k 7 ' 



A^ cosh {(3n^/2) - cosh (f2^[/3/2 - t]) where 



sinh {(3il-~j 



(98) 



fiT)=dDiT) 



dpJT) 



(106) 
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We notice that the function of eq. (|105() does not depend 
only on the relative distance r — r but also from the 
distance of electron from its localization position in the 
crystal. Then the ca. (|32|l becomes 



■j'dr^D.ir) 






pi{r)pi{r',T) 
(Pi(^)>T 



(107) 



We assume f — (electron in its lattice point) and then 
we obtain the variational radial induced charge density 



9Tir) = ^(2.)"-^ f'drD^ir) ^'^'t'l (108) 



2e ' ■' Jo ' "'■'(27r^2(^))3/2 
by eq. H34|l we obtain the variational polaron radius 



can approximate the integral in ea. (|109|l as 

Do{0) f "dT 



drDoir) 

(27r^2(^)) 



^/^ ° Jo {2TrP{T)f^^ 

m 1 / 



2Trh', 



1 - erf, 



drDoir) 



' 2<^(x) 



^ 



(2^(m2))3/2- 

Collecting these results we get eq. (^SJ. 




/3 






drDoir) 



Rp,T= |D^ l^drD,iT)f{T) 



1/2 



(109) 



High density limit 

The characteristic length ^^(r) defined in ea. pO()|) 
is expressed in term of r-dependent positional fluctua- 
tions dD{T), eg. 1)92(1 . which is an integral of a function 
du^ k [t) weighted by the DOS p{Ld) of the Wigner lattice, 
eq. (l98|I .To have an estimate of this integral we replace the 
integration by inserting an average frequency in the func- 
tion d^^ ^. We choose iop/ ,/e^ because it is the typical 
frequency of the electronic fluctuation in the crystal for 
the high density regime eq. ((47|l . Moreover we consider 
the low temperature limit (fc^T ^ ^pj \JEco)- Then 
from ea. H106|l we get the following estimate for 1?{t) 



e{r) 



miopj y 



(110) 



The characteristic time scale of electronic diffusion in 
imaginary time is t^x = iyip j ^£aa) ■ The rising-time is 
Tel = ^IC^^p/ \/£oo)- Therefore we have approximately 



-(r) ~ — r {T<^Tei) 
m 



e{r) 



2m- 



(r > Tel) 



Now in the variational polaron radius Rp^r of eo. H109|l 
another time scale appears Tph = w^o but at high density 
Tph » Tel- Now we can separate the lowest time scale 
Tel contribution in the imaginary time integral so that we 
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